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Abstract 

We investigate propagation of perturbations of equilibrium states 
for a wide class of ID interacting particle systems. The class of sys- 
tems considered incorporates zero range, if-exclusion, mysanthropic, 
'bricklayers' models, and much more. We do not assume attractivity of 
the interactions. We apply Yau's relative entropy method rather than 
coupling arguments. 

The result is partial extension of T. Seppalainen's recent paper 
[firf . For < (3 < 1/5 fixed, we prove that, rescaling microscopic 
space and time by N, respectively N 1+ @, the macroscopic evolution 
of perturbations of microscopic order of the equilibrium states is 
governed by Burgers' equation. The same statement should hold for 
0</?<l/2asin Seppalainen's cited paper, but our method does not 
seem to work for (5 > 1/5. 

1 Introduction 



In the recent paper T. Seppalainen proves that in the so-called totally 
asymmetric stick process (equivalent to Hammersley's process as seen from 
a traveling second class particle), small perturbations of microscopic order 
of equilibrium states, macroscopically propagate according to Burgers' 
equation, if hydrodynamic limit is taken where space and time are rescaled 
by N, respectively N 1+/3 . This result is valid for any < (5 < 1/2 fixed 
and goes even beyond the appearence of shocks in the solution of Burgers' 
equation. Seppalainen's proof relies on the combinatorial peculiarities of 



Hammersley's model and on coupling arguments. It is conjectured in 11] 
that the result should be valid in much wider context, actually Burgers' equa- 
tion should govern propagation of disturbances of equilibria (in this scaling 
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regime) for essentially all interacting particle systems with one conserved 
observable, which under Eulerian scaling lead to a nonlinear 1-conservation 
law. Seppalainen's cited result and also our present paper conceptually is 
closely linked to the work of R. Esposito, R. Marra and H-T. Yau, where 
this kind of intermediate scaling was first applied for the simple exclusion 
model in d = 3. 

In the present paper we partially extend Seppalainen's result. We prove a 
very similar result universally holding for a wide class of interacting particle 
systems. Our proof is structurally robust, it does not rely on any combina- 
torial properties of the models considered. We apply Yau's relative entropy 
method rather than coupling arguments. We pay, of course, a price for this 
generality: (1) applying the relative entropy method, our results stay valid 
only up to the emergence of shocks in the Burgers' solution and (2) we can 
prove our theorem only for f3 € (0, 1/5) instead of the ideal (3 € (0, 1/2). 

Technically speaking, the proof is a careful application of the relative 
entropy method. However, we should emphasize that there is some new 
idea in the 'one-block replacement' step, where the standard large deviation 
argument is replaced by a central limit estimate — and a stronger result is 
gotten. See Lemma ^ and its proof. Also: since in our scaling regime we 
have to consider mesoscopic blocks of size N 2 @ rather than large microscopic 
blocks, in the one block estimate so-called non-gradient arguments (e.g. 
spectral gap estimates) are involved. 

The paper is organized as follows. In section || we present the mod- 
els considered and some preliminary computations (infinitesimal generators, 
equilibria, reversed processes, eulerian hydrodynamic limits, formal pertur- 
bations). In section |3| the main result is precisely formulated in terms of 
relative entropies. Section |i| contains the proof. This is broken up in several 
subsections, according to what we consider a logical structure. 

2 Preliminaries 

2.1 The models 

2.1.1 Notation, state spece 

Throughout this paper we denote by the discrete tori Z/iVZ, N G N, 
and by T the continuous torus R/Z. 
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Let z min , z max G ZU{-oo, 00} with z min < z max , and S := [z min , z max ]nZ. 
The state space of the interacting particles system considered is 



n N -.= s jN . 

Configurations will be denoted 

z := (zj) jeT N G n N , 

2.1.2 Rate functions, infinitesimal generator and examples 

Following ||], Q and § we require that the rate function c : S x S 1 — > [0, 00) 
satisfy the following conditions: 

(A) For any x,y G S 1 

c(^ m i n , y) — — c(x, 2^ max ), 

Note, that this condition is restrictive only if either —00 < z m { n or 
z max < +00. It guarantees that, with probability 1, the local 'spins' Zj 
stay confined within the bounds [z m in, z max ] . In order to avoid degen- 
eracies we also assume that for x G (z m i n , z max ] and y G [-z m irn -^max) 

c(x,y)>0. (1) 

(B) For any x,y, z G S 

c(x, y) + c(y, z) + c{z, x) = c(y, x) + c(z, y) + c(x, z). 

(C) For any x, y, z G S \ {z min } 

c(x, y - l)c(y, z - l)c(z, x - 1) = c(y, x - l)c(z, y - l)c(x, z-l). 

This condition is equivalent to requiring that there exist a function 
r : S — > (0, 00), with r(z m i n ) = 0, such that for any x,y G S \ {z m i n } 

c(x, y — 1) r(x) 
c{y,x-l) r(y)' 

If — 00 < z m \ a or z max < +00, we formally extend r to Z as r{x) = 
for x < z m i n , and r(x) = 00 for x > z ma _ x . 
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Remarks: (1) The monotonicity condition c(x, y + 1) < c(x,y) < c(x + l,y) 
would imply attractivity of the processes defined below. We do not require 
this property of the rate functions. Our arguments do not rely on coupling 
ideas. 

(2) In the case of unbounded z-variable, max{|z m i n |, |.z max |} = oo, we shall 
also impose some growth condition on the rate function c(x,y). See condi- 
tion (D) below. 

The elementary movements of our Markov process are: (zj,Zj + i) — > 
(zj — 1, Zjjf.% + 1) with rate c(zj, More formally, we define Qj : £1 N — ► 

The infinitesimal generator of the process defined on the torus is 

L N f(z)= Y, c(z„z J+1 )(/(e^)-/U)). 

jeT N 

Clearly, due to the nondegeneracy condition (pi) , the only conserved quantity 
of the process is ^ • Zj . 

Remark on notation: Consequently, we shall denote by z = {zj)- eTN 
an element of the state space £l N and by £(s) the Markov process on fl N 
with infinitesimal generator L N . 

There are three essentially different classes of examples. 

(1) Bounded occupation number. The only example with z m \ n = and 
-Zmax = 1 is the completely asymmetric simple exclusion model. For any 
K > one can easily check that there exists a finite-parameter family 
of models with z m [ n = and z max = K satisfying conditions A to C. 
These are usually called generalized K -exclusion models. 

(2) Occupation number bounded from below. There exists an infinite-para- 
meter family of models with z m i n = and z max = +oo. In particular, 
with 

c(x,y) = r{x) = l {x>0} r(x), 
we get the zero range models. 
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(3) Unbounded signed occupation number. From the infinite-parameter fam- 
ily of possible models with z m \ n = — oo and z max = +00 we point out 
the following: let r : Z —* (0, 00) satisfy 

r(z)r( — z + 1) = 1. 

Define 

c(x, y) = r(x) + r(—y) 

Following |IJ, Q we call these models bricklayers models. 

If the occupation number is not bounded (i.e. the state space is not 
compact) we need some additional conditions on the growth of the rates. In 
order to avoid lengthy technical computations we only consider two special 
cases: the zero range model and the bricklayers model, defined in examples 
(2) and (3). For these models we need the following extra conditions: 

(D) (Growth condition for zero range and bricklayers models) 

(i) sup \r(x + 1) — r(x)\ < a\ < 00. 

(ii) There exists xq £ N and 02 > such that r(x) — r{y) > 02 for all 
x > y + XQ. (That means that for x S N r{x) is essentially linear.) 

These conditions will guarantee the existence of dynamics and cf. || the 
uniform spectral gap estimate stated in Lemma [||. 

2.2 Equilibrium states and reversed process 
2.2.1 Stationary measures 

From the growth condition D it follows that 

00 n 00 n 

z = n r{ >- k + i ) + 1 + n r ^) _i < °°- 

n=l k=l n=l k=l 

We define the following probability measure on S 



tt(x) := < 



Z- 1 Y[ rik)- 1 if x>0, 

k=l 

—x 

Z^Ylri-k + l) if £<0. 



k=l 
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For 9 G M let 

F{9) :=logJV*7r(*) 

and 

# min := mf{6 : F(0) < 00} # max := sup{0 : F(9) < 00} 

For 6* G (0min, ^max) we define the probability measures 

vr e (z) :=7r(z)exp{#z-F(0)} 

on S 1 . Expectation, variance and covariance with respect to the measure irg 
will be denoted by E#(- • • ), Var^- • • ) and Covg(- ■ ■ ), respectively. 

According to Q, pH , @, conditions A to C guarantee that for any 
G (^mirn^max) the product measure 

■= n 7r& - 

is stationary for the Markov process. However, due to the conservation of 
Y2j z ji on the finite tori T N these measures are not ergodic. It is a standard 
matter to check that the measures conditioned on the value of ^ • Zj , 

U) : = ^6 U\ ^2 z i = fc )> k €Zn [Nz min , Nz max \, 
j 

are ergodic. We shall refer to ttq , respectively, 7r fc as grand canonical, 
respectively, canonical measures for our model. (The different uses of the 
subscript should not cause any confusion.) 

2.2.2 The reversed process 

The elementary movements of the reversed process are (zj-\, zj) — » (zj— 1 + 
1, Zj — 1) with rate c(zj, Zj-i). 
Define 0* : Q N -> fi", 

= ^ - + ^j-i- 
The reversed generator on the torus : 

L N *f(z)= 53 c(z,,z,_ 1 )(/(9*z)-/(z)). 

Note, that the reveresed process is the same for any ttq , 9 G (# m in, #max), or 
vr^ , G Z n [iV 

^-rnin j A^Z max J • 
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2.2.3 Some expectations 

We denote 

v(9) :=B e (z) = ^zir e (z)=F'(e). 

Elementary computations show 

v \0) = F"(6) = Var e (z) > 0, 

thus (0 m i n ,9 max ) E) i ^ v(ff) G (zmim -Zmax) is invertible. With some abuse 
of notation denote the inverse function by 9{v). 
Further notation: we shall denote 

<5>j := c(z j+1 ,Zj), 

$(v) := E ( v )($j) = ^ ir e{v) {x)ir e{v) (y)c(x,y). 

x,yeS 

Clearly, if — oo < 

■^min ^max ^ then $(t?) is bounded. On the other 
hand, for the zero range models and bricklayers' models with rate function 
r satisfying condition (D), straightforward estimates show that 

< C\v\ 

and also that <l>j has finite exponential moment with respect to any grand 
canonical measure. 

Remark on notation of finite-base cylinder functions: If \E r : 

S m — > R, then we shall denote VPj := ^f(zj, . . . , Zj+ m -i). The indices j £ T N 
are always meant periodically, mod N. Expectation of *$>j with respect to 
the grand canonical measure tt^Lv is denoted 

*(u) :=E 0(v) (%) = ^ tt^^i)...^)^)^^!,...,^). 

zi,...,z m g5 

2.3 Hydrodynamic limits 

2.3.1 Eulerian scaling and its formal perturbation 

For the local density v(t, x) of the conserved quantity ^ ■ Zj, under Eulerian 
scaling, by applying Yau's relative entropy method (see or chapter 6 of 
||, or section 8 of ||), one gets the pde: 

8 t v + d x $(v) = 0. (2) 
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2.3.2 Perturbation of the Euler equation 

Throughout the rest of this paper vq G (^minj Zmax) will be fixed and the 
shorthand notation 

a :=$(v Q ), b :=&(v ), c := $> ) (3) 

will be used. Note that bo is the characteristic speed for the hyperbolic pde 
(Q), corresponding to vq. Furthermore, it is assumed that Co 7^ 0. 

We now consider a small perturbation of the trivial constant solution 
v(t, x) = vq of (§). We fix p > and insert in (||) 

v ie \t,x) := v + e /3 n(e 1+/3 t, e(x - b t)). 

Letting e — > 0, formally the inviscid Burgers' equation is gotten for «: 

ftn+|^(u 2 )=0. (4) 

3 The main result 

3.1 Further notation and terminology 

Let vo € (^min, -^max) be fixed and ao, bo and Co as defined in (||), Co 7^ is 
assumed. We also denote #0 := 9( v o)- 

Furthermore, let u(t, x), t E [0, T], x G T, be smooth solution of Burgers' 
equation (Q). We shall use as absolute reference measure the stationary 
measure 

We define 

6 N (t,x):= N p (9 (vo + N~ u(t , a; - A^&ot)) - #0) 

i.e. 0(u o + N-Pu(t, x - A^M)) = O + AT^tf" (t, x). 

The partial derivatives of 9 N (t,x) are easily computed: 

0^ (t, x) := d x 6 N (t, x) = 0> o + N- p v(t, x - N^bot)) d x u(t, x - N^b t) 

0? (t, x) := d t 6 N {t, x) = -6% (t, x) x (c u(t, x - N\t) + N%) (5) 

In the computation of dtO N we use the fact that u is smooth solution of (|j). 
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The time dependent reference measure (not to be confused with the ab- 
solute reference measure!) is 

v t : = II ^e +N-^e N (t,j/N) = II K6(v +N-Pu{t,j/N-N(>bot))- ( 6 ) 

jei N jeT N 

The true distribution of our process on , at macroscopic time t, i.e. at 
microscopic time N l+ @t is 

mT :=^e^{N^tL N }. (7) 

The Radon-Nikodym derivatives of these last two probability measures on 
Q N , with respect to the absolute reference measure 71^ , are 

f t N U) := ^(z) (8) 
= f[ exp { Zj N-^9 N (t,j/N) - F(9 + N^9 N (t,j/N)) + F(8 )} 

K(z) := ^ r (z)=ex P {N^tL N *}h^(z) 
3.2 What is to be proved? 

We want to prove that if [1q is close to Vq , in the sense of the relative 
entropy H{n$ \ v$) being small, then nf stays close to in the same 
sense, uniformly for t 6 [0,T]. 

How close? Given two smooth profiles U{ : T — > R, i = 1, 2, let 



N 



]^[ n 6(vo+N-0uiti/N))> *— 1,2. 



Then, an easy computation shows that the relative entropy | ) is 

I ) = X] H ( 7T e(v +N-l 3 U2(j/N)) I ^(vo+iV-^uiO'/JV))) 

= N^O'o [ ( U2 - Ul ){u 2 -^(u 2 + Ul ))dx + 0{N 1 -^), 
Jt * 

where 9' := 9'(vq) and Fq := F"{9q). This suggests that one should prove 
ff w (t):= J ff(^| V f)= (JV 1 -^), (9) 
uniformly for t £ [0, T]. 
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3.3 Main result 



Consider a generalized mysanthrope model with rate function satisfying con- 
ditions A-D. Let vq G (Zmin, 2 max ) be fixed so that cq defined in @ is nonzero. 
Let u : [0, T] x T — ► M be a smooth solution of the inviscid Burgers' equation 
(Q). Further on, let uf , respectively, /if be the time dependent reference 
measure, respectively, the true distribution of the mysanthrope process, de- 
fined in M), respectively, (R). 

Our main result is the following 

Theorem. Let (3 G (0, 1/5) be fixed. Under the stated conditions, if 

hOj% \n N ) = o(n 1 - 2 ?) 

and ftQ) holds for t = 0, than ^) will hold uniformly for t G [0, T] . 

Remark: The statement should hold for (3 < 1/2, but, with our method, 
seemingly only /3 < 1/5 can be treated. 

Prom this theorem, by applying the entropy inequality the next corollary 
follows: 

Corollary. Under the conditions of the Theorem, for any smooth test func- 
tion (p : T — ► K 

N ~ 1+P E ^{{j-N l+ %t)/N) (0(iV 1+/3 t) -v ) f <p(x)u(t,x)dx, 
as N — > oo. 

4 Proof 

Our strategy is to get a Gromwall type estimate. We shall prove 

H N (t) - H N (0) <C f H N (s)ds + Err^ (t). (10) 
Jo 

It is assumed that H N (0) = o(N 1 ~ 2/3 ) and the error estimate Err^ (t) = 
o(iV 1_2 ^) is the main point. 

Important remark on further notation: In the remaining part of the 
paper, without loss of generality, we assume 

Vq = 0, O = 0, a = 0. 

This means that from now on z, v, 0, and stand for z — vq, v — vq, 8 — 8q, 
— ao and $ — ao 
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4.1 Estimating d t H N (t) 

In order to prove an inequality like dlCf ) we need to estimate dfH N (t). Using 
the well known inequality 

fLlogfKLf 

which holds for every / > 0, straightforward computations lead to 

d t H N (t)<N^ L^dtf- ?*dtf. (11) 
Jn N It Jn N Jt 



(See chapter 6 of || or the paper |12[ for details.) 

Further remarks on notation: In subsections |4j] and i.1 t € [0, T] will 

be fixed. In order to avoid heavy notations, in these subsections we do not 
denote explicitly dependence on t. In particular we shall use the following 
shorthand notations 

e N (x):=e N (t,x), 6 N x {x):=e N x {t,x), %{x):=B?{t,x), 
U N ( X ) : =u(t,x-N\t) 

Discrete gradient of functions g : T — > R will be denoted 

V N g(x) :=N(g(x + l/N)-g(x)). 

4.1.1 Computation of L N *f t N /f t N 
After straightforward calculations we have 

AT* fN 
I 



jeJ N 

-N- 1 -? 0^j/N)(^-HN^u N (j/N))) 

j£T N 

-N- 1 -? £ 9^(j/NMN-V(j/N)) 

j£T N 

+ Yl (e~ N ~ X ~ B{yN eN){ ~ j/N) - e'^ 1 ^* ^W)®] 

+ AiN-^e^j/N))^ 
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where in the last line the shorthand notation A(x) := e x — 1 + x is used. 
The main term is the first sum on the right hand side. We introduce 

*i := fy-boz, 

:= ^e(v){^j) = Hv)-b v 

and write in the main term 

Thus, eventually we get 

±-J*-dtf = (12) 
Jn N ft 



£ C(i/iV) / fa - N-W (j/N)) dp? 



+ Errf (t) + Err^(i) +Err^(t), 
where the error terms are 
Errf (i) = -J2 0^(j/NMN-^u N (j/N)), (13) 

Err^(i) = iV i+^^( e -^- 1 -' 3 (V JV ^)(i/iV)_ e -iV-^^f y/^) /" $^f, (14) 
ieT iv 

Err- (t) = iV 1+ ^ £ A(lT 1 -f^ (j/N)) [ ^ dtf . (15) 

jeT iV ■>« 

4.1.2 Computation of d t ft I ft 



Now we turn our attention to the second term on the right side of (|ll|). 
From (||) and @ we get: 

= £ iV-^f (j/iV) - AT-^ W (i/iV)) 

- - ^ C (i/iV) (coiV-^^ (j/N) + 6o) - N-?u N (j/N)) 

j£T N 
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In the last sum we write 



and note that the second term is a small error. 
Eventually we get: 

^fW(l) = (16) 
n N ft 

E 6 N X (j/N)V(N~?u N (j/N)) [ [ Zj - N-?u N (j/N)) di4 (z) 
-bo E 6 * ^/N) [ {z 3 - N~?u N (j/N)) dtf + Err? (t) 

where 

Err£(t) = - £ eZUM^iN-WtiM-coN-WV/N)) (17) 

x / {zj-N-WV/Nfidtfiz). 



Note that, when inserting in (|1 1]), the second sums on the right hand side 
of (II4) and (|FJ) cancel out. 

4.1.3 Blocks 

Throughout the paper the one-block size I will be chosen, depending on the 
system size N, so that asymptotically 

/ » N 2 ?. 

We introduce the block averages 

l-l l-i 

i=0 i=0 



The main terms (i.e. the first sums on the right hand side) in (12), respec- 
tively, in (|l6|) become 



- E C 07*0 / w (*} - *(N-?u N (j/N))) dtf + Err? 1 ' (t), (18) 

• gT iV 
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respectively, 



E 8 N x {3/NW{N^u N {j/N)) / $-N-f>u N (j/N)) dtf (z) (19) 

jeT iv 

-n N.l , , n 

+ Err 6 ' (i). 

After rearrangement of sums the error terms Err^'' (i), respectively, Err^'' (t) 
are written as 

Err™ (t) = -£ (r 1 E " (0' " 0/#) " C Cj/^V)) (20) 

jgTJV i=0 



2-1 



Errf (t) = E ( r 1 E # ((? " i)/N)$'(N- V ((j - i)/N)) (21) 



- C (j/N)*'(N-Pu" (j/N))\ / ^ d/i? (*)■ 

4.1.4 Sumup and estimate of the error terms (so far) 

Summing up, from (|TT|) , (|l2|), (ll~6|), ([i~8|) and (|l9|), so far we have got: 

< - E °* (H N ) [ N (*5 - ViN-WU/N)) (22) 

+ §'(JV-/»u w (j/AQ) (4 - AT V (j/iV)) ) dtf (z) 

+ Errf (t)+Err^ (t) + En£(t) 

+ E<(t)+Err^(t)+Err6' ; (t) 

with the error terms given in (|l3|), (14), (p~5|) , (|l7|), (|20|), (|2l|), respectively. 
For the estimate of the these terms we use the following lemma: 

Lemma 1. Let ^ : TL m — > ]R 6e a finite cylinder function and denote := 
^(zj, • • • , Zj +m -i). Assume that, for | — y | < 70; E 7r (exp{7 x I'}) < oo. T/ien 
i/iere exists a constant C < oo depending only on m and 70 , such that for 
any ip N : T N -> K, 

E ^0') / *j < C max |^(j)l (^^ + WEttW), 
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uniformly for t £ [0, T}. 



Proof. We may assume that max JgT iv \ipN(j)\ = 1 and E 7r *I'(£) = 0. We set 
7i := 70-^"^ < 7o then with the entropy inequality: 



gT iV 

< -F( / xf|vr Ar ) + -logE,exp{ 7l £ 



Applying the Holder inequality to the second term, and using that and 
axe independent if \j — k\ > m we have 

I Y, M) [ Vjdtf <-H(^\tt n ) + — ^ A( 7l m^(i)), 



where we use the notation A( 7 ) := log E„- exp{7^(C)}. 

Now, A(0) = A'(0) = 0, thus we have the asymptotics A( 7 ) = 0(j 2 ) 
for ||7|| <C 1. Since max- eT Jv IV'aKj)! = 1 an d 7i = 0(iV _/3 ) there exists 
a positive constant Ci such that A( 7 imV>jv(.7')) < Ci 7 f for every j G T^. 
There also exists a constant C2 with H^nf \ir N ) < C^iV 1 " 2 ^. From these 
the lemma follows with C = C2/70 + Ci^m. □ 

By Lemma |l| and the smoothness of u(t, x) we readily get: 

0(N l -W), 
O(N-0), 
0(N l - 4 P), 

0(N- (3 l), 
0(N- 2f3 l). 

4.2 One block replacement 

On the right hand side of ( p2|) we replace the block average ^?j(z) by its 
'local equilibrium value': ^f(z l j). We denote 

R(x, y) : = *(z) - $(y) - &(y)(x - y) (23) 



Errf (t) = 
Err^(t) = 
Enf (t) = 
Err£(t) = 
Err^ (t) = 
En? 1 ' (t) = 
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Then: 



d t H N (t) < -J2 d *ti/ N )[ R(z l yN-Pu N (j/N))d^{i 
. Jn N 



r N,l 



+M N \t) + 0{N l ^y N-H), 

< sup \e N x {3/N)\ Y, I \R(^N^u N (j/N))\d^{z) 
o<t<T f—t Jn N 



r N,l 



+M N ' l (t) + 0(N 1 - 4/3 V N-Pl), (24) 
where 

M m (t) :=-J2 9 * [ N («} " *C#) drf (*). (25) 



ft i\/r N ' 1 , 



The estimate of J M ' (s)gLs is done in the next subsection, by the so-called 
'one block estimate'. 

We estimate now the first term on the right hand side of (24). Assume 
N = Ml. By the entropy inequality 

+ -log( I ex P { 7 £ |i?(4,iV-V(i/iV))|}^f (. 

We estimate the integral in the second term on the right hand side of (]2< 
using again the Holder inequality: 



exp{ 7 Y \R(z l p N-V(j/N))\}dv?(z) 

gT iV 

/ M-l 

ex P {7^ E \R(4l+i,N-Pu N m + i)/N))\}dv? (z) 

i=l k=0 
r M-l 

J] / exp{/ 7 \R(4l+i,N-V ((kl + i)/N))\} dv? 

i=l JU k=0 

J] / exp{Z 7 |i?(4,iV-V(i/iV))|}^ffe)) 1/ ' (27) 



In the last setp we use the fact that for any fixed i £ [1, i] the block averages 
CL+i' ^ = 0, 1, . . . , M — 1, are independent under the measure . From 
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fl23[ ) it is easy to see that the function 

x ^ R(x + ArV (j/N), N-Pu" (j/N)) (28) 

is asymptotically quadratic if \x\ <C 1. If the variables z% G S are bounded 
than (p8]) is automatically bounded. If S is unbounded, but condition D 



holds, than (28) is asymptotically linearly bounded for |x| S> 1. Thus we 
may use Lemma || stated below, and eventually from (|26|), (|27|) we get for 
7o sufficiently small and I > I/70: 

Consequently, using this bound in ( |24| ) we find 

S^ff* (t) < Cff" (t) + M*'' (t) + 0(A^- 3/3 V AT-^Z V iVr 1 ) , (29) 
holding uniformly for t G [0, T]. 

Lemma 2. Xei ("1,(^2,... fre i. i d. random variables with zero mean. 
Assume 

A(A) := logE (e A< -) < 00. (30) 



Let the smooth function G : K — > 6e quadratically, respectively, linearly 
bounded for \x\ <C 1, respectively, \x\ 3> 1, i.e., < Ci(|x| A (x 2 /2)) ; 

mf/i some finite constant C\. Then there exist constants 70 > and C < 00, 
such that for any < 7 < 70 and I > 1 /70 

Eexp{ 7 ZG((Ci + --- + C/)A)} <C - ( 31 ) 

Remarks: (1) It is worth comparing the statement and proof of Lemma 
^ with the corresponding places in previous works applying the one-block 
replacement, see, e.g., Proposition 1.6. in Part 6. of B. There usually a 



weaker statement (o(l) instead of 0(1) on the right hand side of (31)) is 
gotten by use of more sophisticated tools (large deviation principle instead 
of central limit estimate). Actually, we do need the sharper 0(1) bound. 
(2) The statement is easily extended: imposing more restrictive conditions 
on A(A), the growth condition on G(x) can be relaxed. E.g., assuming 
A(A) = 0(A 2 ) for |A| 3> 1, we may take G(x) quadratically (rather than 
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linearly) bounded at |x| 1. 
(3) Actually, 

lim Eexp { 7 ZG((Ci + • • • + CO/0 } = i 1 ~ 7A"(0)G"(0)) 

l— >oo 



-1/2 



But, since we need only the bound ( |31| ) and not the exact value of the limit, 
we leave the proof of this as a funny exercise for the reader. 

Proof. First we prove the statement with the more restrictive assumption 
A(A) < C2A 2 /2. Assume 7 < (C1C2) 1 and let £ be a standard Gaussian 
random variable, independent of the variables Q. We denote by < ••• > 
expectation with respect to the variable £. Then we have the following 
chain of (in)equalities: 

Eexp{ 7 ZG((Ci + ••• + &) A)} < Eexp{C l7 (Ci + --- + 0) 2 /(20} 



B(exp {^C^/l (Ci + ■■■ + Cl) £}) 

(Eexp {7^7^(6 + ■ ■■+&) e}; 



= <exp{ZA(v^i7Ae)}> 
< <exp{C 2 C l7 £ 2 /2)}> 

= (i- 7 c 1 c 2 )- 1/2 . 

Now we consider the general case. Choose a so large, that for any i£l 

G(x) < lncosh(ax). 

One can do this due to the bounds imposed on G. Let £1, £2, • • • be i.i.d ran- 
dom variables which are also independent of the Q-s and have the common 
distribution P(£j = ±a) = 1/2. We shall denote by < • • • > expectation 
with respect to the random variables £j. We choose Ao, C3 so that for 
|A| < Ao the quadratic bound A(A) < C3A 2 /2 holds and fix 7 < Ao/a. Then 
we have: 

Eexp { 7 ZG((Ci + • • • + Ci)/0 } < cosh (« (Ci + • • • + Ci)/0 ™ ^ 

< E< exp { (6 + • • • + £ M ) (Ci + • • • + Ci)/l} > 
= (Eexp {(6 + • • • + e M )(Ci + • • • + co/0> 

= (exp{/A((6 + ---+£ M )A)}> 

< (exp{C 3 (6 + --- + £ M ) 2 /(2/)}}- 
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Now, since lncosh(aa;) < a 2 x 2 /2, we can apply to the random variables 
the argument of the first part of this proof, with C2 = a 2 and C\ = C3, to 
get 

Eexp{ 7 /G((Ci + --- + O)/0} < {l-jC 3 a 2 y 1/2 . 

□ 



4.3 The one block estimate 

The objective of this section is to provide an estimate for J* M N ' 1 (s) ds, 



where M ' (s) is given in fl25[). 
4.3.1 Cutoff 

We cut off large values of the block averages. In case of compact state space, 
i.e. —00 < z m j n < z nmx < 00 this step is completely omitted. Clearly, 

M N "(t)<A N K '(t) + B N K \t), (32) 

where the terms on the right side are defined as 

A K (t) := 

£ eZ(t,j/N) I W-*(^))i { |*j|vai*ji<jf}^ is), 



B N K \t) , 



sup \6%(t,j/N)\ V / |*5--*(4)|l{^.|va|^.|>^} d ^^)' 

jeT iV jGT iv 

where a > is a fixed constant which will only depend on the rate function. 
For the estimate of Bj£ (t) we need the following lemma (applied with m = 
lor2 only): 

Lemma 3. Let A : Z m — > R 6e a finite cylinder variable. Then there exists 
a map K \— > e{K), such that limx^oo = and 

E L \^-n-W{\^>K } d^{z)<e{K)N^. 
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Proof. The entropy inequality yields: 

<i(F( M fK-)+logE^exp{ 7 £ |¥}-*(<J)|1 { |AJ|>*}}) 

We note that the j th and k th terms are independent in the last sum if 
\j — k\ > I + m — 1. By the Holder inequality, for I > m, we have 

log E^ exp { 7 ^ W-*(Cj)|l{|A$l>*}} 

< JVZ-MogE^ exp{2Z 7 |^ - $(C')|1 { | A «|>a:}}- 
Next we use Cauchy-Schwarz inequality: 
E^exp{2/ 7 |^-$(C i )|l { | Ai | > ^ } } 

<1 + E W * (l {|Ai|>K} exp{2Z 7 |^-%')|}) 

< l + jp^flA'l >K)} 1/2 {E^exp{2/ 7 |^-$(C')|}} V2 . 

From standard large deviation arguments it follows that there exists a func- 
tion [0, oo) 3 7 i— > A(7) G [0, oo) (finite for any finite 7 !), such that 

B nN exp {2l 7 \¥ - $(C')|} < exp (ZA( 7 )}. 

On the other hand, using again a Holder bound and a standard large devi- 
ation estimate, for large I we have 

P w n (\A l \ > K) < mexp{ - U{K)/{2m)}, 

where x ^ I(x) is the rate function 

I(x) := sup (Ax — logE^jv exp {AA}) . 

A 

We define 

7 (K) := sup{ 7 : A( 7 ) < I{K)/{2m)} A K. 

Since lim^oo I(x) = oo, we also have limx-><x> 7 (if ) = oo. Now, putting 
together all our estimates, we get 

£ / |«5-*(4)|i {|A ,| >J0 ^u) 

- 7 T^y( F ^ kJV ) +ivrl(1 + ^ ) )- 
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Noting that H(f$ \tt h ) = 0(N 1 - 2 ^) and I > CN 2 ?, the lemma follows with 
e{K) = C-f(K)- 1 . □ 

It is easy to see, that the functions Aj = zj and Aj = tyj satisfy the 
conditions of the Lemma ||, thus it follows that there exists a map K — > e(K) 
with limii'^oo e(K) = and 

B N K l {t) < e^N 1 - 2 ^ . (33) 

4.3.2 General tools 

We collect in this paragraph the general, model independent facts used in 
the one-block estimate. 

Let C(s) be a Markov process on the countable state space f2, with ergodic 
stationary measure tt. Denote by L and L* the infinitesimal generator and its 
adjoint, acting on L 2 (Q, tt). We denote by D(f) the Dirichlet form associated 
with the generator L and stationary measure 7r: 

D(f) :=- [ fLfd7r = - [ fL*fdir 

The spectral gap of the infinitesimal generator L is p~ l defined by 

P = p(L)-= sup €(0,ooJ. 
/eL 2 (n,7r) ^AJJ 

Actually, this means that (L + L*) /2, the symmetric part of L, has a gap of 
size in its spectrum, immediately to the left of the eigenvalue 0. 
If V : — > R is a bounded measurable function we denote 

ct(L + !/(•)):= sup {spectrum of (L + L*)/2 + V(-)}. 

The following statement is the variational characterization of the 'top 
of the specrtrum' of a self-adjoint operator over a Hilbert space. It can be 
found in any introductory textbook on functional analysis. 

Fact 1. For a (L -\- V (■)) the following variational formula holds : 

a(L + V{-)) = sup(J V(-)hdTT- D(Vh)^, (34) 

where the supremum is taken over all probability densities with respect to the 
stationary measure tt. 
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The second fact is a perturbative estimate of a (L + eV(-)). It can be 
found, e.g., as Theorem 1.1 in Appendix 3 of ||. 

Fact 2. // V : Q — > M has zero mean, i.e., f^V dir = 0, then, for every 
B<(2\\VH ooP (L))- 1 

w(L + eV ^ i T^W^m Vai ' <n (35) 

The third general fact to be used is a direct consequence of the Feynman- 
Kac formula and straighforward euclidean (inner product) manipulations. 
Its proof can be found, e.g., in j/J or as Lemma 7.2 in Appendix 1 of § . 

Fact 3. Assume now that V : M + x f2 — > M. is a bounded function. The 
following bound holds 

E^expj^ V(s,((s)) ds} < exp j J a(L + V(s,-)) ds}, (36) 

where now denotes expectation over the Markov chain trajectories started 
from the stationary initial measure n. 

4.3.3 Notations 

We shall use the notation lx , respectively, fj! for a generic probability 
measure on £l N , respectively, ! . We shall denote by h N (z), respectively, 
h l (z) their Radon-Nikodym derivatives with respect to the absolute reference 
measures 71"^, respectively, ir l . Further on fx ' ' 3 will denote the [j, . . . ,j + l — 
1] marginal of fi N and \i ' := iV -1 ^2j e jN M ' ' 3 the average Z-dimensional 
marginal of fx N . Correspondingly, h " 3 (z), respectively, h ' (z) will denote 
the Radon-Nikodym derivatives of [i " 3 , respectively, Li ' , with respect to 
the absolute reference measure it 1 
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For k £ Z fixed we denote: 

l 



n* ■■= 


i=l 






m k := 








i 

w k : = 


A*'( n k)> 








i=i 


*) = 


-1 ^ 


/4(l) 


l 

■■= ^ (= Yl Zi = 

i=i 


k)- 






. it /4U) 


= 1 





Denote by D N , D ( respectively D l k the following Dirichlet forms 



D N (f) 


1 N f 

i=l ^ 


-/U)) 2 ^U) 




i N r 




D\f) 


1 i_1 /■ 

:= oE / c («i>«i+i)(/( e i«)- 






i=2 ^ 


-/U)) 2 ^(z) 


D l k (f) 


i=i ^ Q fe 


■ /(l)) 2 




= ^/ i c(^,^ 1 )(/(e* 1 )- 

i=2 ^fc 





In the definition of D N periodic, in that of D l and D l k free boundary condi- 
tions are understood. 

It is easy to check that for any probability measure fi l on Q l 

^W)=E^U\/>4)- (37) 

Further on, using convexity of the Dirichlet form one can readily prove that 
D N (VhF) > j D^y/hF™). (38) 
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4.3.4 Applying F-K formula 

We return now to the concrete computations. Before the estimate of 
J * Aft 1 (s)ds we need some more notation (we do not denote explicitly de- 
pendence on the cutoff): 

V"' l (z) := (*} " *(^))l{|*»|Va|*«|<K-}» 

V\z) :=V?' l {z), 

V^ l {t,z) :=e»(N-Wh,j/N)V;- l (z), 

V N - l (t,z):= ^ V 3 N '\t,z). 
jeT N 

We denote by £ N (t) the Markov process on Q N with infinitisimal generator 
L N and by E^at , respectively, E^at the path measure of this process starting 
with initial distribution /j,q , respectively, 71"^ . 

By the definitions and the entropy inequality we have 



ft 1 /*iv "Z 

J Q 4F(°)*> = j^^(j o V N >\s,f(s))ds] 

<_i TF ( ff (^|^) + logE 7rN exp{ J lV N ' l ( S ,( N (s)) ds}). 

We apply the Feynman-Kac bound ( |36l) and the variational formula (34) to 
the second term on the right hand side of the last inequality: 

log E^ exp { J -fV N ' 1 (s,C (*)) ds] (39) 

< / a(L N + 7 V N ' l (s,-))ds 
Jo 

sup ( / 7 V N ' 1 (s, -)h N dn N - D N (Vh? )) ds. 
h N y Jn N ' 
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Using (^) we bound the integrand in the last expression 

supf f 7 V N ' l (s,-)h N - 1 dn l -D N (VW)) (40) 



= sup( £ [ iVf'Mh™ <rf -D»{y/iFj) 

hN jeT N Jn ' 

<supV ( [ 7 V"'\s,-)h N ' l ' j d-K 1 -IdUV^ 

< - V sup ( [ IjV*' 1 (s, -)ti dit 1 - D l (Vh l 
1 .tit h« Wn< 



jgT JV /i 



Next we use (0) and again the variational formula (| 

sup ( f ljV*'\s, -)h l d-K 1 - D l {VU )) (41) 

supj^u;[sup( / l-yV 3 N ' l (s,-)ti k dn l k - D l k (Jh l k 
w l k h' 



sup ^2 w k a { L k + hVj' (s, •)) 

S u P ^^(z 7 C(^i/^)EU^) 

+ a(4 + / 7 C( s ,i/iv)(^-EUu i 



In the first step we used (|37|), The second step is a straightforward identity. 
In the third step we have used again (|34| ) and we introduced the notation 
L' k for the infinitesimal generator of the process restricted to £l l k . Finally, in 
the last step we use the notation introduced at the beginning of the present 
paragraph. 

4.3.5 Spectral estimates 

The rest of the proof of the one block estimate relies on the following three 
steps: (1) a straightforward estimate of E l k (V l ) and Var^(V ! ); (2) a lower 
bound of order ~ l~ 2 on the spectral gap of L l k , valid uniformly in k G Z; (3) 
combining these two and the perturbational bound (35), an upper bound on 
(x(. . . ) appearing in the last expression. 
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Lemma 4. There exist constant C(K) < oo for every K > Kq, such that 
for any I and k the following bounds hold: 



\K(V')\ <C(K)l 



-1 



Var'JV 1 ) <C(K)l 



-i 



(42) 



Proof. For \k\ > Kl there is nothing to prove, so let \k\ < Kl. Restricted 

i 

k 



on n l 



V l =¥- *(k/l) - - *(fc/0) >K y 



Consequently, 



\K{v l )\ <2 



+ E i ,f|^ i -E'^'ll 



\>K} 



)■ 



By the equivalence of ensembles (see e.g. Appendix 2. of H and also 



-1 



The second term can be estimated with the Cauchy-Schwarz inequality and 
with large deviation techniques (noting that because of the growth condi- 
tions on the rates we can choose such a > that aT l K > |E^' | uniformly 
for \k\ < Kl) and it can be easily shown to be smaller order then the first 
one. Var^(y') may be estimated with similar methods. □ 



Lemma 5. There exists a constant C < oo, independent of I and k, such 

A) 

;2 W 



that for any f G L 2 (Q' k ,7r' k 



Vax l k {f)<Cl 2 D l k [f). 



(43) 



Proof. For the details of the proof of this gap-estimate we refer to || , || , 
||. For models with bounded z-variable, — oo < z m \ n < ,z max < oo, we note 
that 

c(x,y) > ar(x)l {x>ZmlBilKw} . 

with some positive constant a. Thus, it is sufficient to prove the gap estimate 
for the reversible process with rates c(x,y) := f(x)Hr x>Zinliit y <Zxugx \, which 
has the same ergodic stationary measures 7r l k as our original process. For this 
latter process the induction steps of || or Appendix 3 of |j] apply without 
any essential modification. 
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In m the statement is proved for zero range model with rate function 
satisfying condition (D). Minor formal (but not essential) modifications of 
that argument yield the result for the bricklayers' models with rate functions 
satisfying condition (D). □ 

Remark: Actually we could consider a wider class of models with un- 
bounded spin space, by imposing 

inf c(x, y) > a r(x) 
y 

with some positive constant a and r{x) obeying condition (D). 

We remark that there exists a constant C depending only on the solu- 
tion u(t,x) of the Burgers' equation (Q), and another constant C(K) which 
depends also on the cutoff level K, such that 

sup \0Z(t,j/N)\ < C, (44) 

0<t<T 
j€T N 

H^-E'fcV'H^ < C(K) (45) 

Now, combining @, @, @ and @, we get the following 

upper bound, which holds for every sufficiently small 7: 

*(L> k + lj9»(s,j/N)(V>-K(V 1 ))) < ^fj^p 

Setting 

7 := 7o r 3 with 7o <min{l, (2C 2 {K))~ 1 } 

we have 

a(L l k + W»(s,j/N)(V l -E l k (V 1 ))) < C(K) 7 2 l' 3 . 

Collecting all the estimates and going backwards through ([ill), ( fig) , (|39|), 
we find eventually 

log exp { J jV N ' 1 (s, C N (*)) ds} < C(K) 10 N 2+ H- 4 



and 



f A N K l (s) ds < C(K)(N- 3(3 l 3 + Nr 1 ) (46) 
Jo 
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Consequently, from (^), and (p2[) , with any fixed K < oo we have 

f M N '' (s) < e(i^)0(iV 1 - 2/3 ) + C(if )(iY~ 3/ ¥ + iVr 1 ) (47) 
Jo 

where C(K) is a finite constant which may increase to infinity as K — > oo, 
and e(-ftT) — > as K — ► oo. 

4.4 End of proof 



We put together fl29|) and ( (47|) to get, for any if < cxd fixed (with a C not 
depending on K) 

H N (t) < H N (0) + C f H N (s)ds + e{K)0(N 1 - 213 ) 

Jo 

+0(N 1 ~ 3/3 V N-^l V iVr 1 V N- 3 ' 3 l 3 ) . 

If 

0</3< i 
5 

then we can choose 

N 2 ? « / « M 1+ ^ 3 

which ensures 

o(n 1 -'^ V N~ P l V Nr 1 V iv~ 3/ ¥) = o(N l ~ 2 P). 
Thus for every if < oo 

H N (t) < iif^ (0) + C I' H N (s)ds + e(if )^- 2/3 + (iV 1 - 2 ^) , 

JO 

where limi<r_>oc e(K ) = 0, and from Gromwall indeed (0) follows, uniformly 
for i G [0,T]. 
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